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Experimentally the charmonium hyperfine splitting is Mj^^— M^^ = 117 MeV and current lattice 
results are generally below this value. The difference could be due to the effects of the dis- 
connected flavor singlet diagrams which have not been included in these calculations. Previous 
attempts to determine the disconnected flavor singlet corrections have led just to rough estimates 
in the range of ±20 MeV. We present preliminary results for these corrections calculated on fine 
{a ^ 0.09 fm) Asqtad 2+1 flavor lattices provided by the MILC Collaboration. 
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1. Introduction 

The lattice calculation of the hyperfine splitting in charmonium is currently still showing a 
discrepancy with the experimental value of 117 MeV. The discrepancy is large (30-40%) in the 
quenched case [|l|, 0]. It is reduced to around 10% in dynamical studies with improved actions [||, 
but its origin is still uncertain. A possible explanation could be that even the current state-of-the- 
art lattice action formulations do not reproduce the quark dynamics within the charmonium states 
sufficiently accurately for the purpose of this calculation. Another possibility which could account 
for the discrepancy (or at least part of it) is the neglected contribution of the disconnected diagrams 
in the lattice computations. These diagrams could contribute to the masses of both the vector J/y/ 
and the pseudoscalar rjc, and thus affect the hyperfine splitting M//^ ~^ric' Perturbatively, the 
contribution of these diagrams in charmonium is expected to be small due to the OZI suppression, 
especially for the vector state. However, non-perturbative effects, such as the C/a(1) anomaly 
and glueball mixing, might enhance it enough so that it becomes a non-negligible fraction of the 
hyperfine splitting. 

In this work, we present our preliminary results of the effort to test this possibility and deter- 
mine the contributions of the disconnected diagrams. Previous calculations [^, ^ using two-flavor 
gauge ensembles very roughly estimate the contribution to be within ±20 MeV. They both con- 
firm that there are significant difficulties in obtaining a signal for the disconnected diagrams due 
to noise, especially for heavy quarks. Our calculations are performed on an ensemble of 505, 2+1 
flavor MILC lattices generated with the improved Asqtad action [^. For this ensemble the 
ratio of the masses of the light and heavy quarks is ruud/^s = 0.1. The charm quarks are simulated 
using the Fermilab interpretation of the clover action with /c^ = 0.127 tuned to the physical charm 
quark mass. The disconnected diagrams are calculated stochasticly with 12 spin and color diluted 
sources. 

We employ several methods for improving our calculation in comparison with previous works. 
First, our lattice volume is 40^ x 96, which is much larger than in [^, ^] and our lattice scale, 
a ^ 0.09 fm, is finer. In addition to these improvements we use the unbiased subtraction technique 
[ llO| ] in the stochastic estimators. The success of this technique depends on the fast convergence 
of the hopping parameter expansion used in the subtraction. Considering that Kc is small for the 
charm quark, we use the terms of the expansion only up to third order in which reduces the 
standard deviation of the disconnected correlator by about a factor of 4. The last and most impor- 
tant improvement is that we study the point-to-point (PP) disconnected correlators instead of the 
traditional time- slice-to-time- slice (TS) ones. As a result, the standard deviation of the signal is 
reduced about three orders of magnitude. 

2. Lattice formulation 

The disconnected part of the flavor singlet TS correlator is calculated as: 

D{t) = cr (L(0)L*(0) , where L{t) = Tr(rM-^) (2.1) 

and the trace is over the Dirac, color and space indices. For the vector we have F = 7^, cp = 1 and 
for the pseudoscalar F = 75, cr = — 1. To determine the effect of the disconnected diagram on the 
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masses of the charmonium states, previous works ^ ^ explore the ratio of the disconnected to the 
connected TS correlators at zero momentum: 

^-^_l-4/,K-/)^_l (22) 

c{t)-c{t) '-a/ ^^-^^ 

In the above F{t) = D{t) +C(/) is the full propagator corresponding to a state with "full" mass 
Mf. The mass rUc extracted from the connected propagator is the mass which is usually studied 
in lattice simulations. The constants Af and are the full and connected propagator amplitudes, 
respectively. Considering that the available lattices are quenched with respect to the charm quark, 
an appropriate fitting form for the ratio data at zero momentum would be 

= ,23) 

The differences nic — mf should be calculated for both the vector and the pseudoscalar and then 
subtracted from the mass of each in order to determine the effects of the disconnected diagrams on 
the hyperfine splitting. 

In this work however, we calculate the PP disconnected correlators, and the above analysis 
has to be modified accordingly. The disconnected PP propagator as a function of the Euclidean 
distance on the lattice r is defined as: 



^W^TT (M^)^*(^0>, (2.4) 



where the sum is over all pairs of lattice points at this distance, A^^ is the number of these pairs 
and there is trace only over spins and colors in L. The connected PP correlator C(r) is defined in a 
similar manner. It is known that in the continuum limit the asymptotic behavior of C(r) is: 

C(r)^A^, (2.5) 

and consequently the behavior of D{r) in this case can be deduced as follows: 

d e~'^''' 

D{r)^---^C{r)^B^. (2.6) 

Thus their ratio is 
D(r) B 

where we interpret the amplitude ratio as: 
B 

-^Mc-mf. (2.8) 

The behavior of D{r) and C(r) on the lattice will be affected by discretization artifacts; still, for 
large r and fine lattices we assume that the above behavior is a good approximation. 
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3. Results for D{r) 

From the previous studies it is known that the charmonium correlator signal disappears very 
quickly around t = 2 — 3. We work with the PP disconnected propagator, since this way we ben- 
efit from both the additional data at non-integer distances and the much improved statistics. The 
correlator D{r) has from one to three orders of magnitude smaller relative errors than the TS dis- 
connected propagator in the region where we have a signal. Figure 1 illustrates this statement by 
comparing D{r) and D{t) for rjc for two different ranges of r and In the right panel of Fig. 1, 
the comparison is done on a shorter range in order to emphasize the fact that we do have a clear 
signal for D{r) in the range where the D{t) signal is completely obscured by the noise. The result 




Figure 1: Comparison of D{t) and D{r) for rfc for two different ranges of t and r. 



that the D{r) signal is so much better than the one for D{t) can be explained by the fact that in the 
calculation of D{t) there are a great number of contributions from points, which although not far 
from each other in the t direction, are far in the 4d Euclidean space. For the disconnected correlator, 
the noise increases strongly with the distance and such points contribute nothing to D{t) but noise. 
This problem is avoided when working with D{r) instead. 

4. The D{r) /C{r) ratio for t]^ 

Figure 2 shows our data for the ratio D{r)/C{r). We make the following observations: 

• The roughness of the data at short distance is a sign of rotational symmetry breaking and 
requires a nonasymptotic lattice fitting function. 

• The ratio D{r)/C{r) inherits the sign flip of D{r). The sign flip is expected, since from the CP 
symmetry of the operator, the sign of D{r) at r = must be opposite the sign of any single- 
pole contributions, and they dominate at large r. We believe this is the first observation of 
light state contributions to an operator constructed from charm quarks. 

• In a separate analysis we find that excited states in the TS correlator C{t) contribute more 
than 20% for r < 5. The same is expected for D{r). Indeed, if we assume a single heavy 
mode at short distance, to make the nonasymptotic form of D{r)/C{r) resemble the data, we 
need a mass much heavier than the T]^, as indicated in the right panel of Fig. |^ The light 
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states in D{r) are very likely dominant for r > 4, but we estimate that the charmonium signal 
is still about 30% of the whole signal at r = 5. 

• The asymptotic expression Eq. (2.7) is not likely to be applicable anywhere. At large distance 
D(r) is dominated by lighter hadronic modes (the T] and T]' mesons, glueballs, etc.). Thus to 
extract the ratio of ground state amplitudes we must resort to more elaborate models. 

• In the right panel of Fig. § we present two more exploratory approximations for D{r)/C{r) 
that include one heavy and one light mode, just to show that the behavior of the ratio data 
can be approximated in this manner. 

5. Extracting the ground state t]^ signal from D{r) 

Following the discussion in the previous section, at present we are led to explore a simplified 
model for the intermediate range r G [5, 11] that includes a ground and excited T]^ state and a light 
state. The light state represents a number of possible states. We work separately with D{r) and 
C(r), extract the amplitudes of the T]^ contributions, and finally use Eq (2.8) to get an approximate 
mass shift. 

For the connected correlator we have only a connected- r/^ ground state and one excited state. 
For the disconnected correlator our ansatz includes a light state as well. In momentum space it 
reads 



2 



^ (5.1) 



p^ + mj 

where m/ is a light state mass, m* is an excited connected rjc mass, and f{p^) and g{p^) are real 
functions. We need the residue of the double pole at + = to determine the mass shift. In 
coordinate space the asymptotic form is then 

f2 fl fl 

The terms correspond, respectively, (Fig. |3|) to double pole contributions from the ground and ex- 
cited charmonium states, mixed ground-excited state contributions, and the light state contribution. 
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Figure 3: Charmonium disconnected diagrams included in the fitting form Eq. (5.2). The ground state 
diagram is to the left, the excited is in the middle and the "mixed" one is to the right. 



For this preliminary analysis, we take the amplitudes (B) of the ground and excited charmonium 
states to be equal. We do this in order to reduce the number of free parameters. Since the excited 
state contribution is just a small correction in the chosen fit range, this restriction on the amplitude 
is not of great importance. The coefficient c is determined from the mass terms in the ansatz. We 
fix c ^ 7, but find that our results are not very sensitive to varying it over the range c G [2, 14]. 

We adjust only B and L. We fix the masses as follows: The light mass is taken to be the 
central value m/ = 0.43(1) from a single-exponential fit (as in Eq. (2.5)) in the range r G [7, 12]. 
The connected T]^ and T]* masses, nic = 1.1598(7) and m* = 1.51(5), are known from fits to the 
connected TS propagator C{t). 

Before fitting we smoothed the data by averaging the signal in small bins in r, thus reducing 
the effects of the lattice artifacts and improving the continuum approximation. The resulting fit to 
D{r) is shown in Eig. H. The fit has X^/^f = 24/19 . Our fit result favors the following range of 
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Figure 4: The disconnected PP correlator D{r) and fit to it in the range of distances re [5,11]. 



values for the ratio: 



— ^mc-ruf e [-4,-1] MeV. 



(5.3) 



Thus the disconnected diagram contribution slightly increases the r]c mass. 
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6. Summary and conclusions 

We have calculated the disconnected pseudoscalar propagator with dramatically reduced noise 
using unbiased subtraction methods and point-to-point correlator data. We fit the data to a sim- 
plified model and extracted the amplitudes of the ground state in the disconnected and connected 
propagators. We find that the disconnected diagram contribution increases the T]^ mass by 1-4 MeV. 

The vector meson J/y/ can be studied in a similar fashion. If the corresponding mass shift 
of the J/y/ is negligible as expected 0], we would then conclude that the hyperfine splitting is 
reduced slightly. 

This result is preliminary, however. Further work is needed to explore the sensitivity of this 
conclusion to the choice of the fitting model and to include effects of rotational symmetry breaking 
properly. Calculations at a smaller lattice spacing are needed to test for heavy-quark and hard-gluon 
discretization errors. A similar study of the //y/ could test the common expectation that its mass 
shift is negligible. 
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